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Abstract. We present hadron-level predictions from the Monte-Carlo generator CASCADE and 
numerical calculations of beauty quark and quarkonium production in the framework of the kj- 
factorization approach for LHC energies. Our predictions are compared with the CMS experimental 
data. 
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^j' Our study is motivated by very recent measurement of open beauty quark and b-']Qi 

production performed by the CMS Collaboration. It was observed that the data tends 
to be higher than the MC@NLO predictions and that the shape of the pseudo-rapidity 
distribution is not well described by MC@NLO. The p^-spectra of b-jets are not well 



Cs| , discribed too [1,2]. 

^T} \ Recently we have demonstrated reasonable agreement between the ^y -factorization 

predictions [3] and the Tevatron data on the Z?-quark, bb di-jets, 5+- and D-meson 

Q . prodiction and also agreement with total set of HERA data for 7/yA-mesons [4]. Based 

on these results, here we give the analysis of the CMS data in the framework of the 
kT-factorization approach. We produce the relevant numerical calculations in two ways: 



We will perform analytical parton-level calculations (which are labeled as LZ). 
The measured cross sections of heavy quark production will be compared also with 



c^ \ the predictions of full hadron level Monte Carlo event generator CASCADE [5]. 

The basic dynamical quantity of the fcj -factorization approach is the unintegrated 
kj-dependent) gluon distribution (UGD) s^{x,kj,ii^) obtained from the analytical or 
numerical solution of the BFKL or CCFM evolution equations. 

The cross section of any physical process is calculated as a convolution of the partonic 
cross section d and the UGD ■{2^(x, fc|,/i^), which depend on both the longitudinal 
momentum fraction x and transverse momentum kj : 

/oo oo oo oo 

=c/g(xi,fcir,/i )£/g{x2,k2T,IJ- )dgg{xi,X2,kij',k2T,---)dxidX2dk^jdk2T- 

The partonic cross section a has to be taken off mass shell (kj- -dependent). 

It also assumes a modification of their polarization density matrix. It has to be taken in 



so called BFKL form: 
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Concerning the UPD in a proton, we used two different sets. First of them is the KMR 
one [6]. The KMR approach represents an approximate treatment of the parton evolution 
mainly based on the DGLAP equation and incorpotating the BFKL effects at the last 
step of the parton ladder only, in the form of the properly defined Sudakov formfactors 
Tg(k^,/i^) and r^(k^,/i^), including logarithmic loop coorections. 

Second UGD is the CCFM one. The CCFM evolution equation have been solved 
numerically using Monte-Carlo method [7]. In this case UGD are determined by a 
convolution of the non-perturbative starting distribution £/q{x) and CCFM evolution 
denoted by £/{x, k^ , /i^) . The parameters of £/q (x) were determined in the fit to F2 data. 

It is known that the hard partonic subprocess of gluon-gluon fusion gg — )• QQ ampli- 
tude is described by three Feynman's diagrams. 

In the numerical calculations in the case CCFM UGD we have used two different sets, 
namely AO and 50. The difference between these sets is connected with the different 
values of soft cut and width of the intrinsic kj- distribution. 

For KMR we have used as the standard GRV 94 (LO) [8] as MSTW [9] (in LZ 
calculations) and MRST 99 [10] (in CASCADE) sets. The UGD dependes on the 
renormalization and factorization scales jAr and jif. We set /i| = Mq + (pjj, +P2y)/2, 

lip = f -f Qy, where Qj is the transverse momentum of the initial off-shell gluon pair, 
ruc^ 1.4 ±0.1 GeV, m^ = 4.75 ±0.25 GeV. We use the LO formula for the coupling 
Ofj(juJ) with nf = 4 active quark flavors at Aqcd = 200 MeV, such that ttsiM^) = 0. 1232. 



We begin the discussion by presenting our results for the muons originating from 
the semileptonic decays of Z^-quarks (Fig. 1). To produce muons from Z?-quarks, we 
first convert Z?-quarks into 5-mesons using the Peterson fragmentation function with 
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FIGURE 1. The pseudo-rapidity distributions of muons arising from the semileptonic decays of beauty 
quarks. The first column shows the LZ numerical results while the second one depicts the CASCADE 
predictions. Solid, dashed and dash-dotted, dotted histograms correspond to the results obtained with the 
CCFM AO, BO and KMR UPD. The experimental data are from CMS. 



default value e^ = 0.006 [11] and then simulate their semileptonic decay according to the 
standard electroweak theory taking into account the decays Z? — )■ /i as well as the cascade 
decay b —^ c —^ jj.. In CASCADE calculations also Peterson fragmentation function is 
used but with full PYTHIA fragmentation. In Fig. 2 we show our discription of the b-jet 
distributions at LHC. 

In the case quarkonium production we used Color-Singlet (CS) photon-gluon or 
gluon-gluon fusion in the framework of the ^j- -factorization approach. In Fig. 3 we 
show our results for 7/1// production at HERA and in Fig. 4 we present comparison 
of our results with LHC data on the 7/1// production in framework of the MC generator 
CASCADE. 

In summary we have analysed the first data on the beauty and 7/i// production in 
pp collisions at LHC taken by the CMS collaboration. Our study is based on a semi- 
analytical parton level calculations and a full hadron level MC generator CASCADE. 
The overall description of the data is reasonable. In most of the distributons it is simi- 
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FIGURE 2. The double differential cross sections da/dydpr of inclusive b-jet production as a function 
of Pt in different y regions calculated at y/s = 7 TeV (LZ predictions). Notation of all histograms is the 
same as in Fig. 1. The experimental data are from CMS. 
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FIGURE 3. Differential cross sections J/y/ mesons at HERA. The solid, dashed and dash-dotted 
histograms correspond to the results obtained using the CCFM AO, BO and KMR gluon densities. The 
upper and lower dashed histograms represent the scale variations. 



lar to MC@NLO except in some particular distributions where the fc^ -factorization ap- 
proach does describe the data better, like in Z?-jet.7/i// production in the fc^ -factorization 
approach with CS model comes much closer to the data than the coUinear calculations. 
The reason is the off-shell ME, which includes even higher order contributions than the 
NLO coUinear calculations. 
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FIGURE 4. Differential cross sections J/y/ mesons at LHC in CASCADE. The experimental data are 
from CMS. 



